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THEORETICAL ANALYSIS OF OXYGEN DIFFUSION AT STARTUP IN AN 
ALKALI METAL HEAT PI PE WITH GETTERED ALLOY WALLS 
by Leonard K. Tower 
Lewis Research Center 

SUMMARY 

The diffusion of oxygen into, or out of, a getter ed alloy ejqposed to oxygenated alkali 
liquid jnetal coolant, a situation arising in some high temperature heat transfer systems, 
was mialyzed. The relation between the diffusion process and the thermochemistry of 
oxygen in the aUoy and in the alkali metal was developed by making several simplifying 
assumptions. The treatment is therefore theoretical in nature. However, a practical 
example pertaining to the startup of a heat pipe with walls of T-111, a tantalum alloy, 
and lithium working fluid illustrates the use of the figures contained in the analysis. 


INTRODUCTION 

Refractory metal alloys have been proposed for use at high temperatures in struc- 
tures for space power production and other applications. When alkali metal coolants are 
present, oxygen in the system components must be minimized or segregated if corrosion 
is to be avoided. One means of segregating oxygen in refractory metal alloys has been 
the use of a gettering metal. To getter oxygen, hafnium has been introduced into 
tantalum-8 tungsten (wt. %) alloys to form T-111 (tantalum-8 tungsten-2 hafnium). Be- 
cause of the stability of hafnivun dioxide (Hf 02 ), most of the oxygen present in, or enter- 
ing, the alloy will be sequestered as Hf 02 until all the hafnium is consumed. 

In some applications involving gettered alloys, the oxygen levels to which the alloy is 
exposed may vary sharply. For instance, in a new heat pipe with alkali metal working 
fluid, the refractory metal wicking, arteries and walls contain oxygen introduced in all 
the fabrication processes. At startiq), the liquid metal purified by distillation from the 
evaporator condenses on the condenser walls and is entrained in the wicking and arteries. 
The chemical potential of oxygen in these parts is very likely to be higher than that of the 
oxygen in the liquid metal , and oxygen will be leached into the liquid and transported to 


the evaporator. Because the liquid metal vapor can contain little oxygen, the concentra- 
tion in the liqviid metal of the evaporator may rise. As this concentration rises, the ox- 
ygen can move into the evaporator walls. If the concentration of oxygen in the evaporator 
wall becomes sufficiently- high over a period of time, corrosive failure of the wall is a 
probability (ref. 1). 

Refractory alloys containing tantalum and niobium can dissolve considerable amounts 
of oxygen (ref. 2) and diffusion at elevated temperatures is rapid (ref. 3). The presence 
of a gettering element as an oxygen sink in alloys of tantalum and niobium can be expected 
to modify both the rate of diffusion with time, and the oxygen capacity of the alloy. 

In this report a simple analysis of the oxygen diffusion process into, or out of, a 
gettered alloy exposed to oxygenated liquid alkali metal is presented. The analysis re- 
lates the thermochemistry of oxygen in the liquid metal and in the container metal, to the 
diffusion into or out of a semi -infinite alloy. The analysis is intended to clarify some of 
the features of the diffusion processes process which may be encountered in cooling sys- 
tems using gettered refractory alloys and alkali liquid metals. For this purpose, an ex- 
ample is presented for oxygenated lithium in contact with T-111, a gettered tantalum al- 
loy chosen because it is being considered for use in heat pipes. The analysis is part of a 
continuing study of gettering systems for use in heat pipes. 


ANALYSIS 

Figure 1, from reference 4, shows a computed distribution of oxygen between lithium 
and tantalum under conditions of equilibrium. At any given point on figure 1, oxygen in 
the tantalum at atom fraction XQ(Ta) same chemical potential as oxygen in the 

dilithium oxide (LigO) at mole fraction x in lithium. Experimental determinations of the 
soliibility of oxygen in lithium and in tantalum, extended analytically, were used in ther- 
mochemical calculations to construct figure 1. Figure 1 must therefore be regarded as 
largely theoretical. While the distribution shown has not been substantiated by exper- 
imental data, it will be used here for purposes of illustration. 

The atomic fraction of oxygen in tantalum against the mole fraction of Li20 in solu- 
tion with lithium is shown for several temperatures. Also shown in figure 1 is a crossed - 
hatched area representing the zone of stability for HfOg computed as though metallic haf- 
nium were in a phase separate from tantalum. The data for HfOg were replotted from 
reference 5. Actually, the hafnium in dilute solution in an alloy such as T-111 will have 
a chemical potential somewhat different from the pure hafnium metal (ref. 6). The effect 
of this change in the chemical potential of hafnium on the location of the hatched area is 
no doubt of less consequence than some of the other assumptions entering into the con- 
struction of figure 1 (ref. 4). 
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In the simple model of a diffusion event considered here, the alloy is assumed to be 
of semi-infinite extent in the u direction (fig. 2(a)). This assumption is made in order 
to simplify the mathematics. The initial concentration of oxygen throughout the alloy is 
c^, taken here in units of density, or grams per cxibic centimeter. At time t = 0 con- 
ditions external to the alloy are suddenly changed, such that the concentration of oxygen 
required for equilibrium in the alloy is c^.' Exposure to oxygen -bearing lithium could 
bring about such a change in conditions. If c^ exceeds the oxygen concentration Cg at 
which the gettering metal oxide is stable, the oxide is formed by the diffusion of pxygen 
at a rate governed by the diffusion constant D. 

The assumption is made that c^ is below the concentration required for oxidation of 
the main alloying constituents (tungsten and tantalum in T-111). Any formation and pre- 
cipitation of the oxide of the gettering metal is assumed to have negligible effect on'the 
diffusion coefficient D. When the concentration Cg is reached at any point in the alloy, 
the getter is assumed to be oxidized at a sufficiently rapid rate that chemical equilibrium 
is maintained locally. 

The method of attacking the present problem was adapted from techniques used by 
Crank (ref. 7) for mass diffusion in which the diffusion constant is discontinuous with 
concentration. This method is similar to that used by Danckwerts (ref. 8) for the general 
case of heat transfer with a moving boundary, as in the freezing of ice. 

The time -dependent differential equation for diffusion in one dimension is 


8c 

8t 



( 1 ) 


Referring to figure 2, one solution pertains to the region Cg < c < c^, designated by 
subscript 1; 


Ci = c„ + A erf — H_ (2) 

2yl^ 

The other pertains to the region < c < c , designated by subscript 2; 

o 

C 2 = + B erfc — (3) 

2-y^ 


where c is the oxygen concentration, u is the distance coordinate , t is time, erf z 
designates the error function, and erfc z = 1 - erf z is the complementary error 
function. 
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Figure 2(a) presents a schematic view of conditions during the precipitation of get- 
tering metal oxide. A lamina of thickness Au is shown at the concentration Cg, within 
which precipitation of getter oxide is occurring. The flux of oxygen into the lamina is 
-D(3cj/9u) and the flux out is -D(3c2/3u), both in grams per square centimeter per sec- 
ond or comparable units. If the reaction between getter X and oxygen is 

mX+nO-*X^Oj^ (4) 

Q 

and the concentration of getter in the alloy is c^^^ (g/cm ) , the oxygen gettering capacity 
of the alloy is 


* — ^ 

q = -Cm- — (g/cm^) (5) 

exclusive of the oxygen dissolved in the major alloy constituents. The negative sign in- 
dicates the consumption of free oxygen, and M denotes molecular weight. 

If q Au grams of free oxygen are consumed in the lamina, the oxygen consmned 
must be represented by the net diffusional influx into the lamina during the time span At, 
or 

( dCf 

D — i-D— i At = qAu (6) 

. 0u 0u/ 

The lamina must move with velocity V = Au/At to encounter fresh getter material. 
Thus equation (6) in the limit gives a condition for the boundary between regions 1 and 2 : 



At the gettering concentration c_ occurring at u=u , both c, in equation (2) and c« 

5 & 1 ^ 

in equation (3) must equal c . Therefore equation (2) becomes 

o 

u 

c„ = c + A erf — i- (8) 


and equation (3) becomes 
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( 9 ) 


u. 


c„ = + B erfc — i— 

g °° I— 

2yjTyt 


Since c is constant, the variable u /2‘>^ must be constant. Defining 

o & 




2^^ 


( 10 ) 


enables equations (8) and (9) to be solved for A and B as follows: 


A = -i 2 

erf K 


( 11 ) 


c - c 

B =-i f 

erfc K 


( 12 ) 


Equations (2), (3), (11), and (12) with a little manipulation can be expressed as 


erf 


u 


°1 ' = 1 . ^ 2\/Dt 


o °° 


erf K 


(13) 


erfc 


u 


^2 - 




2 \^ 


- c_ erfc K 

O ■ oo 


(14) 


where 


a = 


c - c 
_0 g 

c - c 
o <« 


(15) 


and 


c - c 

P = -2 f = l-a 

■''cc 


( 16 ) 
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Eqmtion (13) applies only when Cj > c^, and equation (14) only when Cg < c^. The get 
tering front velocity from equation ( 10) is 



(17) 


Combination of equations (7), (10), (13), (14) and (17) results in an implicit ejq)ression 
for K 


2 

^+y\^Ke'^ =0 (18) 

erf K erfc k 


where 


r = 



(19) 


Essentially, y represents the relative gettering power of the alloy since it is the ratio of 
the oxygen gettering capacity to the uptake of uncombined oxygen required to achieve 
chemical equilibrium. Equation (17) can be solved iteratively for k, with values of /3 
and y assigned. 

-1 2 

Figure 3(a) shows k for several values of /3 against y from 10 to 10 , while fig- 

2 S 

ime 3(b) shows k against y from 10 to 10 . As the relative gettering strength y in- 
creases, K decreases, representing a slowing of the gettering front (eq. (17)). 

With K determined as a function of )3 and y, equations (13) and (14) can be eval- 
uated. Lettii^ 


z =^y/m 
2 


(20) 


equations (13) and (14) give 


c - c 


c - c 
o «> 


oo 1 a erf z n / „ ^ 

— = 1 - 0 < z s K 


erf 


( 21 ) 


c - c 


oo o erfc z 


= /3 


c - c 
O oo 


erf 


K < Z < 00 


( 22 ) 
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Figure 4 shows the concentration function (c - c^)/(c^ - c^) plotted against z for /3 of 
0, 0.2, 0.4, 0.6, and 0. 8 and values of y ranging from 10“^ to 10^. 

As y; representing the gettering strength, increases, the break in the curves of 
(c - c^)/(Cq - c^) becomes more pronounced, as required by equation (7). When /3 = 0 
as in figvire 4(a), no material passes through the gettering front. For this reason the 
case /3 = 0 applies to one -dimensional diffusion in bodies of thickness y vintil the get- 
tering front reaches y. 

The total amount of oxygen contained in the structure at any time, per unit area, is 
given by 


= Ac.du.r 

Jo Ag 


c« du - q u 


(23) 


where, Cj is given by equation (13) and C 2 by equation (14). Integration of equation (23) 
and combination with equations (10), (18), and (19) gives 


bi = 2 + 2y^ (24) 

(Cq - c^) ^f^ erf k 

A plot of the dimensionless gettering parameter G/2'^Dt (c^ - c^) against y for several 
values of 0 is shown in figure 5. 

The foregoing development treated the diffusion of oxygen into an alloy. The same 
analysis and figures can also be applied to the outgassing of oxygen from a material con- 
taining uniformly distributed oxygen in excess of that required to oxidize the getter. In 
such a case c^ > c^ > c^, and q (eq. (5)) is positive since the getter acts as a soimce 
rather than a sink. Figure 2(b) illustrates the situation. Figiires 3, 4, and 5 are again 
applicable, with G/2\^ (c^ - c^) in figure 5 pertaining to loss of oxygen from the body. 


EXAMPLE 

. The following example illustrates an application of the figures that have been pre- 
sented. The example is based on the situation that might occur in the evaporator of a 
new heat -pipe. The refractory metal wick structure and arteries of alkali metal heat 
pipes often contain large amovints of oxygen after fabrication. At the first startup of the 
pipe, a large amount of oxygen from these sources and from the condenser walls is swept 
into the evaporator causing a sudden, large increase in the oxygen concentration of the 
alkali metal in the evaporator section of the pipe. 
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Consider an alloy with a composition within the range designated as T -111 (ref. 9). 

By weight, this alloy contains approximately 89.5 percent tantalum, 8 percent tungsten, 

2. 5 percent hafnium, 0. 003 percent oxygen, 0. 003 percent nitrogen, and 0. OOl percent 
carbon. A semi -infinite section of this alloy is assumed to be heated to 1700 K and si- 
multaneously exposed to solution of lithivim and lithivun oxide, with the mole fraction x 
of LigO in solution equal to 0. 02. The problem will be to find ^ and y and from these 
determine the time for the gettering front to move 1 millimeter and the time for 
(c - c^)/(Cq - c^) to reach 0. 5. The amount of oxygen per square centimeter of struc- 
ture when the front has penetrated 1 millimeter is to be found. The equilibrium distribu- 
tion of oxygen between liquid lithium and tantalum in the alloy is to be taken from 
figure 1. 

The assumption is made that any oxygen in the system must either be dissolved in 
the tantalum, or combined as Hf02 and dissolved in the tantalum. The oxygen concen- 
trations must be below the point at which tungsten oxides or TagOg forms. The possible 
formation of lower oxides of hafnium is disregarded. The possible presence of neg- 
ligible amoimts of nitrides is also ignored. 

The criteria for the distribution of oxygen in the alloy as fabricated are the 
followii^: 

(1) When Cq < c^, no HfOg is present and c^ = Cq. Here C denotes the total 
concentration of an element in the alloy, both combined and uncombined, after 
fabrication. 

(2) When Cg < Cq < 2Cjj^ (Mq/Mjj^), some amoimt of hafniiun is oxidized and 

c^ = Cg. Moreover the effective gettering capability is reduced because some of the haf- 
nium is already oxidized during fabrication. The amount of fabrication oxygen which 
must be taken up by hafnium is Cq - Cg. This reduces the gettering capability of the haf- 
nium, and equation (5), modified to account for the fabrication oxygen, becomes 

-q = 2C„ ^ - (Cq - V (25) 

WlRf 

(3) When Cq > 2Cjjj (Mq/Mjjj), all hafnium is oxidized and no further gettering 
occurs. 

The density p of T-111 is nominally 16.70 grams per cubic centimeter (0. 604 lb/ 
cu in. ) (ref. 9). The total concentration of hafnium for the assumed composition, both 
oxidized and unoxidized, is 


CRf = ^Hf^ = 0. 4 18 g/cu cm 


(26) 
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where w is the weight fraction of constituent in the alloy. The total concentration of ox- 
ygen is 


Cq = Wqp = 5, 02x10“^ g/cu cm (27) 

From figure 1, for 1700 K and a LigO mole fraction x of 0. 02, x^ = 3. 1x10" for 
monatomic oxygen dissolved in tantalum. Converting this to an oxygen concentration in 
the alloy gives 


a R 

^ - g/cu cm 

g M ^ 

Also 2Cjj| (Mq/Mjjj) = 0. 075 gram per cubic centimeter from equation (26). Criterion 2 
is thus satisfied, and c^ = Cg = 4. ixlO"® gram per cubic centimeter. From equa- 
tion (25), -q = 0. 074 gram per cubic centimeter. 

The oxygen concentration c^ represents equilibrium between the oxygen in the tan- 
talum and in the lithium after the hafnium has all been consumed. From figure 1 for 
X = 0.02 and T = 1700 K, = 6X10"®. Then 

6 

“=o = 'Ta =^(Ta) = 
iviTa 

Therefore, c^ - c^ = 3.8x10" gram per cubic centimeter. From equation (16) /3 = 0 
and from equation (19) y = 1. 95x10^. For these values of /3 and y figure 3(b) gives a 
K of 5X10"^. 

In order to compute the time for penetration of the gettering front to a selected depth 
from the foregoing value of y, the diffusion coefficient for oxygen in T- 111 must be 
specified. Reference 3 presents this equation for the coefficient of diffusion of oxygen in 
tantalvim 


D = 0. 0044 exp 


12 800\ 
T ) 


cmVsec 


(28) 


-3 

At 1700 K, D = 8.474x10 sqviare centimeter per hour. This value is accepted here for 

-3 

the T-111 alloy. Since k = 5x10 , the time for the gettering front to reach 1 millimeter 

(eq. (10)) is 1.2x10^ hours. 

Q 

From figure 4(a) a ratio (c - c^)/(c^ - c^) of 0. 5 is reached at z = 2. 5X10" . For 
the diffusion coefficient just determined, the ratio of 0. 5 is reached at 1 millimeter of 
depth after 4. 7x10*^ hours, by equation (20). This value at first glance seems to 
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represent an inordinately slow diffusion process in view of the 1700 K temperature. 
Three factors contribute to this effect. One is the gettering process which causes a 
steep concentration gradient to the left of the front, and a very flat gradient to the right 
(fig. 2). The other factors, the oxygen distribution coefficient from figure 1 and the as- 
sumption of a single diffusion mechanism, will be discussed further in the next section. 

For i3 = 0 and y = 1. 95x10 figure 5 gives a value of 295 for the gettering param- 
eter G/2(Cq - c^)V5t. After the front has penetrated 1 millimeter, the total amoimt of 
dissolved and combined oxygen in the structure is then 1. 26x10"^ gram per square centi- 
meter. In view of the very large value of u almost all of this oxygen is contained as 
HfOg between the liquid -solid interface (u = 0) and the gettering front (u = Ug). 

While the example herein was inspired by a problem involving a lithium -filled heat 
pipe with finite getter ed walls, the figures that were used pertain to a semi -infinite body. 
In this example, a value of /3 = 0 was determined which meant that no diffusing oxygen 
was found on the side of the front away from the oxygen source. In such cases the semi- 
infinite solution applies exactly until the gettering plane passes entirely through the 
material. 


DISCUSSION 

In the foregoing analysis and example, numerous assumptions were necessary in 
order to achieve a simple and understandable solution. Real physical situations can be 
expected to deviate considerably from this ideal situation. 

The distribution of oxygen between lithium and tantalum used in the example was de- 
termined from fig^e 1, taken from reference 4. While figure 1 was based upon exper- 
imental determinations of the solubility of oxygen in lithium and in tantalum, thermo - 
chemical theory was used in its construction. A good experimental determination of the 

distribution coefficient has not been made. However, if the actual ratio was within the 
3 8 

same range, 10” to 10" , as the theoretical value on figure 1, it might defy measure- 
ment. If the actual value of the distribution coefficient were much larger than the one 
determined from figure 1 for the example, the gettering front would progress at a pro- 
portionally higher rate than was calculated in the example. 

The simple theoretical treatment discussed in this report wotild become more com- 
plex if possible ternary oxides of lithium, tantalum, and oxygen are formed. Ternary 
oxides are suspected as part of the corrosion process in contacts between refractory 
metals and alkali liquid metals (refs. 10 and 11). While they have not been observed 
with tantalum and lithium, there is always a possibility that they could play a role in the 
intergranular penetration of oxygen into the refractory metal. This could alter the dis- 
tribution coefficient (fig. 1) and the diffusion coefficient (eq. (28)). 
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Further complexity in this ideal treatment would be introduced by the formation of 
HfgO and HfgO (ref. 12) in addition to the HfOg assumed, Fiirthermore, the solution of 
large amounts of oxygen in hafnium (refs, 12 and 13) could cause additional problems. 
Whether this can occur when the hafnium itself is a minority constituent of a solid solu- 
tion is problematical. The present analysis could possibly be modified to incorporate 
the additional gettering implied by solution of oxygen in hafnium and the formation of ^ 
lower oxides. 

The assumption of local chemical equilibrium at the gettering plane led to equa- 
tion (7). Although in some cases the getterii^ front might move so rapidly as to inval- 
idate such an assumption, in the example just considered, the front moved very slowly 
after a small penetration and local chemical equilibrium was justified. 


CONCLUDING REMARKS 

The diffusion of oxygen into, or out of, a getter ed alloy exposed to oxygenated liquid 
metal coolant was analyzed. Such a situation arises in heat transfer systems having get- 
tered refractory alloy tubing with alkali liquid metal coolant. 

In the analysis, the relation between the diffusion process and the theoretical ther- 
mochemistry of oxygen in the alloy and in the alkali metal was developed. Certain as- 
sumptions such as the semi-infinite extent of the alloy metal wall were necessary to 
achieve a simple solution. For these reasons the treatment is largely theoretical in 
nature. 

The use of the figures developed in the analysis was illustrated by an example for 
oxygenated lithium coming in contact with T-111, a gettered tantalum -based alloy, at 
high temperature. This situation may be encountered in the startup of some heat pipes 
for high-temperature applications, which have been proposed and tested. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 12, 1973, 

503-25. 
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Atomic fraction of oxygen in tantaium, XQ^-g) 



Mole fraction of U2O in solution with lithium, x 


Figure 1 - ~ Equilibrium distribution of oxygen between liquid lithium and solid tantaium at severai temper- 
atures with zone for stable Hf02(c). 



(bl Diffusion out of alloy. 

Figure 2 . - Schematic illustration of oxygen diffusion into and 
out of semi-infinite section of gettered aiioy. 
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(a) Range of y, 10'^ to 10‘ (b) Range of 

Figure 3. - Gettering plane position parameter k slotted against y for several values of p. 
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(a) p = 0. 

Figure 4. - Variation of concentration function (c - c„)/(Cq- C oo) with z ” u/2Vot for severai vaiues of y. 
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(b) p - 0. 2. 

Figure 4. - Continued. 
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(Op-0.4. 

Figure 4. - Continued. 
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(d)P"0.6. 

Figured. - Continued. 
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(e)& = 0.8. 

Figure 4. - Concluded. 
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NASA-Langley, 1973 
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Figure 5. - Gettering parameter G/2>/Dt (c^ - Coo> plotted against y for several values of p. 
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